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io  ei spicy  tho  orcirs cry  notation  of  tensor  analysis,  e„  g„, 

Aj-v  :  tensor  field  of  rank  5,  contravariant  of  rank  2, 

153  corariont  of  rank  3,  all  indices  running  independent¬ 

ly  from  1  to  3. 
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Christoffel  symbols  of  the  second  kind 


al  a2 ’ • *  ak 
**i  ^2  *  *  * 


if  any  two  superscripts  or  any  two  subscripts 
are  equal  or  if  the  &re  nob  the  same  numbers 
as  the  b^; 


if  b-,  ».  b.  may  be  obtained  from  a*  u2  aj.. 

by  an  even  perniutation* 

if  b-y  .  b^  ray  be  obtained  from  a?  a^. 


We  employ  also  tho  special  notations  of  Murnagban  (ref  „  n) ,  which  are 
explained  in  section  IX,  Xn  particular,  Indices  written  to  the  right  denote 
F uteris i  tensor  components  -  while  indices  written  to  the  loft  denote 
Lagrangian  tensor  components.  For  example,  denotqo  the  covuriant 
derivative  with  respect  to  the  Isgrangian  ob&rdimie  la  of  tho  Kulerian 
contmvariant  vector  field  s1 , 


I,  Introduction 

A  aewhouian  continuum  we  shall  define  an  a  finite  or  Infinite  closed  re 
gicr  of  Euclidean  thrse-  dimensional  space,  at  each  point  of  which,  with  the 
possible  exception  of  isolated  sits  of  points  of  dimension  less  than  throe, 
there  exist  the  following  •riaitive  characteristic  functions j 

density  p  j 

stress  tenser  ? 

.internal  energy  per  unit  nans  e; 

specific  entropy  a; 

he at- flow  vector  Hj{ 

end  whose  topological  tronofornation  in  tisse  into  a  finite  number  of  other 
closed  regions  is  governed  by  the  principles  of  conservation  of  rings  rooaen- 
turn,  and  energy  ,, 
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t  it  pirp'S-.s  nre:  v  }  to  state  the  throe  fu-^mcoial  pv5 nek  pi-; s  in 
their  rest  g*jn9'."%l  forsaj  (8)  to  Calculate  the  eowitionr;  gcvot  ’ling  i-srf  .ctly 
elastic  bodi-3.»»  {0}  to  fvmiirts  tee  aquations  g.  ,reriu:?:?  vis.  our>  f laids. 

Our  traatiscnfc  1,-;  not  r jst ’icteri  to  homogeneous  or  isotropic  todies,  nor  arc 
»9  eoccsraad  .itb  the  linear  a] prcoclre.ti.ona  employed  in  the  classical  theories.. 
In  put  treat®  .at  of  perfectly  elastic  solids  V3  follow  ihtfnaghah  (ref.  u), 
whose  notation  rx.  adopt,  venarrlisirg  his  results  so  that  therp  apply  fee 
nonho»oger-30U5  and  non -m:  If  oral y  her  feed  bodies..  Pros  our  point  of  vie®  a  dls- 
tLiguiehing  e  laxcterictie  of  hurnaghan's  work  is  that  ho  defines  an  o.la$t.i .? 
solid  act  by  t.r,f  particular  atracr-strain  relations,  linear  o:  otherwise,  hi' 
by  specifying  t-i&t  the  free  energy  function  must  depend  only  upon  the  gradients 
of  the  deforr-tion  and  that  tfco  temperature  he  constant  during  the  deformation; 
he  thus  ;■  v«?s  r.-atbpmticei  expression  to  cur  intuitive  concept  of  perfect 
Blr.irticity*  fie  due?  leal  linear  Hooke’s  law  then  appears  in  his  theory  as  «; 
first -order  &  atroxi nalieu  2ho  pain  object  of  the  preeaat  paper  is  to  present 
a  similar  for  lulatior.  of  the  theory  of  viacoue  fluids 0  We  aeok  a  imthcai.it i- 
oat  espreasioa  for  the  intuitive  concept  of  a  fluid  and  are  led  to  speciiy 
tb  j  mniKn*  in  t/i  iC  dissipation  function  my  depend  upon  the  velocity  fi  a* 

dio.'its,.  the  pressure  and  the  temperature;  the  classical  linear  Savior* 

3tni.es  I'.w  appears  In  our  theory  aa  a  first- order  approximation,  Oar  gaiwo  e). 
dissipation  function,  if  broken  off  after  the  cubic  terms,  becomes  identical 
dth  the  VAncoua*  portion  of  the  dissipation  function  tor  slip  flow  cs*l«U" 
is  ted  by  durnstt  (quoted  in  raf»  b,  p,.  656)  frost  the  kinetic  theory  of  ganec 


F inegatica 


bet  ‘a  Pc  the  initial  or  L&grttngian  coordinates  of  a  point  In  an  arbi  • 
V’JTj-  ruclidt  an  coord ins  ce  system,  f  nd  x1  the  final  or  tulsrian  coordinate  a 
c.  tic  same  point.,  Thu  motion  of  the  material  coptijmuf  is  expressed  eh  «. 
s>'curnce  of  •*.pe],o$>x.i  brans 'orris  taocs,  s*  ■?  :r  ( !a,  a,;,a,  t), 

-  x.  >%:•?*  t),  •where  t  ic«  the  time-  A  moving  point,  identified 

t  •  the  :xgrR*).J*;iAn  cocrdl.mUva  ra  ••7l.ll  often  bo  called  <?.  wirticle.  The  co< 
p:>:  snt. 5,  of  t'ie  I/*.gvar  gj.t-  a  •awl  .'blevlan  strain  tensors  p7«  aid  ^  ore  df  - 
rerpec hr.'we  ly  c-Sc?  feeds'  as  ' 


{ r&/ 


/  /  «  V'  '  ! 

v  l  a:  Mr-  V  , 


/*t«r  r /•(  f  ...  <  c  /lx"  l 

Sio/  -  ,-y  ••  v-  --7J,  -V'  wi¬ 
lt  O’?.  ho  the  .vpin.-x  trleaent  of  arc  leugrh  at,  «t  specified  initial 

b^  ?r>d  ( .ic )  la  tnat  at  tire  t, 

I!  *'  vo.vp^^nfcs  of  tb  w.(lcolzv  vector  V  &ro  ^iven  \y  tho  tibi 'initios 

v  -  ,  ,;>  / 

V ■  '  r  V1  /  -a  ,4i; 


..  lbv>.  ft  A 


/?  =  (  ^7Z~)  +}  i\^\/Y 

^*<vv  -  u, 

giveS’^fS^JfiSuS  ^/0r  °f  a  function  of  «*  Sulerian  variables  is 

w  ^  •  (?) 


v  ~7F 


'■-  6jC 

•*»  -<■<!  I 'll  ■  ■»  ^ 

<£/ 


exi^^rr  °thGr  ^*«rL  ^  AnLvSt^  Sth'th,  aci!ier,auon» as  ^ 

^S; 

exce^ic  f?d?Ce  "*  ”°b  to  '0e  expected  to  Jcma^r^lld6^^^1  ftquatloE^  which 

”sS2^rtian>-  "ws  ^  -*  *»  *4^T«*tSfss  cr  s:r 

donn^or^^  °f  tha  rate  Sensor  g  y  w  f:i-,en  bj  ^ 

%*£( Xj+K<)- 

Hurmghan  (refc  a,  p.243)  has  shewn  that 

"{>  +  £*  ‘A  )  ,  (8) 

to  13  not  ^  motion  re  inti? . 

euetnrary  to  approximate  the  cemponantu'o?  t£ee*iCfty  J)  St  is 

by  ths  compcnerte  of  tbs  .  It‘^®  £^5&£&i3Sto  to;-, rot  /f.-/-,t 

*““*««  th.  t»  oeS  of  teSS  **  S?*W  imit L- 


£< ;  -  o 


(9)‘*mw* MtSiea  «j'TboS“Jl°„fillt6  lnt°rral  *f  «■>»  th0  eoasttto 

4r^’  —  r 

*t  ~ °  no 
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throughout  tha  interval,  It  is  possible,  howavor,  that  a  bod;,-  may  be  in 
rigid  motion  at  e.  given  instant ,  although  at  instants  before  and  after  it  is 
deformed.  The  cr  iterion  for  instantaneously  rigid  motion  cannot  b©  the 
equation  (10),  since  the  components  of  the  3t,rain»rate  tensor  are  by  equa¬ 
tion  (8)  functions  of  the  components  of  the  strain  tensor,  which  are  computed 
with  respect  to  a  fixed  initial  configuration  which  can  in  no  way  influence 
the  kinematics  of  the  present  instant,  k  necessary  and  sufficient  condition 
for  instantaneously  rigid  motion  is  given  by  the  equations  of  Killing  (ref,  i 
P.244): 

Ci/  -  0  •  (11) 

The  components  of  th©  deformat ion-rat©  tensor  are  in  fact  the  time  rates  of  .• 
change  of  the  components  of  the  strain  tensor  when  th®  referonce  configura¬ 
tion  Mfith  respect  to  whj.ch  the  components  of  the  strain  tensor  are  computed 
is  th©  infinitesimally  jt seeding  one;  that  is,  if  £  y  (tq,t^)  are  the  com¬ 
ponents  of  the  strain  tensor  computed  at  the  time  with  respect  to  the  con¬ 
figuration  st  time  t^,  than 

jr..  -Li*,  ;i2) 


6t 


ZL  i  fry-  — 

Ai-*c 


it  £■</  jfi,  3 

At 


The  law  of  force  of  a  perfectly  elastic  body,  which  responds  only  to  its  de¬ 
formation  from  a  preferred  initial  state,  we  shall  thus  expect  to  employ  tho 
components  of  the  strain  tensor  The  law  of  fore©  of  a  type  of  plastic 

body  which  responds  both  to  a  preferred  initial  state  and  to  th©  rate  at 
which  it  is  being  deformed  from  tiiat  state  we  shall  expect  to  employ  both  the 
<rerr  prevents  of  the  strain  tensor  &  and  the  components  of  the  strain  rate 
tons nr  rf  0  The  law  of  force  for  a  fluid,  which  exhibits  no  response 

whatever  to  any  preferred  state  but  resists  being  instantaneously  deformed , 
to  shell  expect  to  employ  only  the  components  of  the  defornation-ratn  tensor 

Bij£ 

III,  Conservation  of  Mass 

Several  different  mithems  tical  expressione  of  the  principle  of  conserva¬ 
tion  of  mass  are  used  in  the  mechanics  of  conticna  (ref.  a,  pp,  244^246; 
ref.  fc,  pp,  142-1>8);  of  these  we  shall  require  two: 


t>^ocv~< Truism  j 
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*jf  t-  =  o  . 


In  these  expressions  /O  is  the  present  density;  the  density  at  time  t0j 
Ii»  Iji  I3  respectively  the  sum  of  the  one-,  two-,  and  three-rowed  principal 
minors  or  the  strain  matrix  computed  with  respect  to  time  t0: 

21s i (« £; e;  j  r,  =  k% I .  <16> 

IV.  Conservation  of  Momentum 

Let  T1^  be  the  components  of  the  stress  tensor  and  F*  the  components  of 
external  foroe  par  unit  volume.  Then  in  the  absence  of  an  external  moment 
field  the  principle  of  conservation  of  momentum  states  that  (rsf .  a,  pp.244- 
246) 


VJ=T’'-  , 


~7~<Q<  e~  < 

T  .«■*- r 


A‘  . 


Let  vba  any  scalar  function  of  the  Eulerian  coordinates  and  the  tia»s 
Components  Wy  of  a  fluid  stress  tensor  are  given  by  the  definition 

W/  r  W  j  *  7"/  '  {] 

The  equations  (17)  and  (IS)  now  become 

W!>~  W'u .  (: 


V/?  ■+■  ‘T„  4-  F-  -  , 


These  results  are  independent  of  what  particular  function  tt  is  used  in  equa¬ 
tion  (19).  -?e  shall  see  later  that  it  is  usually  convenient  to  let  77"  be  the 

thermodynamic  pressure  Jh  ,  to  be  defined  presently,  although  sometimes  it 
is  the  stress  pressure  ffgiv^n  by  the  definition 


We  postulate  that  the  energy  aeeooiated  with  a  continuum  nay  be 
divided  into  two  parts*  Internal  enfyar  £  per  unit  naes,  and  kinetic  or 
yV2  per  unit  naes.  Potential  energy,  if  It  estate  at  a3lj  we  regard  as 
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energy  of  the  field  of  body  force  against  which  the  particle  does  work,  not 
m  energy  of  the  particle  ltcolf»  Thus  the  totnl  energy  of  a  f  ini  to  volume 
of  the  continuum  is  fcfas  integral  cf  p(&  v-£  \/)  over  the  volume  a  then  the 
principle  of  conservation  of  energy  states  that  the  rods  at  vrlrj.ch  energy  is 
lost  by  an  arbitrary  finite  volume  of  ths  continuin'*  Is  equal  to  the  rate  at 
which  perk  is  dene  against  body  forces  plus  tho  rate  at  which  the  stress  ten¬ 
sor  does  work  upon  the  boundary  surface  plus  tbs  rate  at  which  heat  energy 
flora  out  across  the  boundary  surface : 


'  i  '  i  i 

where  are  the  components  cf  the  heat.;- flow  vector  por  unit  aroa_,  measured 
in  Bochanical  units.  In  problems  of  sisip'le  conduction. 


(24) 


where  T  is  the  temperature  and  k  the  coefficient  of  heat  conductivity,  but 
wo  shall  not  use  this  ausunption ,  lotting  K.  be  sufficiently  general  to  in« 
elude  heat  leases  by  radiation,  etc.  T.x  m  apply  tl«  divergence  theorem  to 
the  surface  integrals  in  forr.ua  (23) »  perform  the  differentiation  on  the 
lofts  and  equate  the  inuegrund  in  the  resulting  voluire  integral  to  zero,  we 
obtain  the  forarjua  of  Aokeret  and  of  do.jal  (ref,  dj  ref,,  e,  p,  5); 


p -ff  -bp  V ~  F  *  +  T  il  %  ^  V  H  *  *  •  (25) 

S  Irani  if ying  this  equation  with  the  aid  of  the  dynamical  equations  (12)  ar;i 
(11* )  -.72  find  that 


Je 


/j>  _ 
t  jy 


?  2  ''  Jr  JV'~  /?’ 

•«?  -V. 


(26) 


This  .a  the  msehanioal  fore  of  the  principle  of  conservation  of  energy.  An 
equivalent  fonaula  has  been  given  by  Jslargulea.,  Lash,  and  Stewart  (ref,  b, 
p.,c4d;  ref,,  1,  p.lSK), 


»7e  now  sketch  the  pstn  ution  construe?, ion  of  a  dbenoiaeiiological  theory 
of  thorixjdynainifls  whit*,  while  loss  clirectly  motivated  physically,  on  prefer 
because  cf  ita  mthcBaticcl  rUiplieity  aod  clarity  to  the  nore  convent  tonal 
expositions  employing  differential  ircra:cnts  ara  aorv  iawod iE.tr If  fAtniliaz- 

primitive  concepts. 


The  thsnnocy mlo  principle  of  conserve tier  cf  energy  ie  an  independent 
postulate  stating  that  the  internal  energy  cf  a  given  {article  say  bo  con¬ 
sidered  at  all  timet  a  function  of  its  density  er.5  a  ien  prirdlvive  scalar 
function  of  the  *a  end  t  which  wo  six!'.!  call  the.  gre s:Lf 1 c  entropy ..  sw 
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(27) 


&(,)  ‘  ,  fy*)  )  —0> 

where  the  parameters  B(i),3/2)  are  given  functions  of  the  lagrangiau 

variables  hx  only.  The. 30  parameters  'to re  introduced  apparently  by  V,  Bjerkner. 
(ref,  f,  pp.  81-33)  to  represent  continuous  inhomogeneity,  such  as  a  distri¬ 
bution  of  salinity  in  water  or  humidity  in  air.  A  discontinuous  inhomogeno- 
ity,  such  as  surface  of  separation  between  water  and  mercury,  la  not  represent¬ 
ed  by  Bjerknes's  parameters,  but  appears  in  the  mathematical  treatment  as  a 
surface  of  discontinuity  across  which  the  differential  equations  are  replaced 
by  suitable  transfer  conditions.  The  equation  (27)  is  called  the  caloric 
equation  of  state:  specifying  the  form  of  the  function  specifies  the  physi¬ 
cal  nature  of  the  continuun.  If  the  parameters  Bn do  not 
actually  occur  in  the  equation  (27),  the  continuum  is  nocoge needed 


The  equation  (27)  is  formulated  on  the  Lagrangian  plan.  To  obtain  a  cor¬ 
responding  Eulerien  equation  of  state,  we  differentiate  the  equation  (27)  k 
times  with  respect  to  tine  (ref,  f,  p,82),  and  then  eliminate  the  parameters 
Bf.-L)  from  the  resulting  kyl  equations,  obtaining  an  end  formula  of  the  typs 


Thus  the  Eulerian  equation  of  state  is  not  a  scalar  equation  but  &  differen¬ 
tial  aquation  of  k— -  order. 

Let  the  specific  volume  be  /: 


'■'=/r 


The  temperature  T  and  the  thermodynamic  pressure 
variables  given  respectively  by  the  definitions 


are  two  new  state 


k.  =  (M  \ 

7 - ‘A 


(29) 


(30) 


(31) 


The  thermodynamic  jjreasurs  ®  does  not  necessarily  have  any  connection  what¬ 
ever  with  the  stress  pressure  p  given  by  the  definition  (22).  We  shall 
henceforth  suppose  that  9  arxi  s  arc  always  functionally  related,  so  that 
the  definition  (30)  is  alway3  meaningful. 

Lot  C  be  a  given  curve  in  the  ©,v,*a,2ft,3a  space,  whose  poiuts  ax# 
parametrised  with  the  parameter  .  If  we  differentiate  the  equation  of 

state  (27)  along  C  atxJ  ejrploy  the  definitions  (30)and  (31)  of  pressure  and 
temp  return  we  find  that 


‘3223 
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NOlk 


A  . 

AS  -  t 

Jp  T 


if  r  %  p  («■  A  (m; 

If  ms  restrict  our  atte  ntion  to  curves  0  along  which  'i  {'Clj^x.  >{  “  6> 
tliat  in,  to  curves  representing  changes  undergone  by  a  given  particle,  we  re¬ 
duce  the  equation  ’  32;  to  the  stapler  fora 


rtf.  \ 


=  =^) 


?  ’  <£$ /'A j  T 

In  other  Korda ,  thu  f  crania 


.jE/Vk  ! 


r<  /  <k  v-  A 

•^7  -  v=r  ?„<.  ?  T" 

A ,  *  o  ^ 


&> 


(33) 


(34) 


Hence 


(35) 


T-  Z-h 

f  s4  v 

is  valid  in  any  con.tinr.usi,  homogeneous  cr  no 
$JL  -  -  h  . 

'tf  * ft 

The  equation  (35)  13  the  usual  postulated  for^  oi“  the  f.iret  law  of  thento- 

dysaid.es,  which  we  here  regard  instead  as  t.  consequence  of  the  postulated 
existence  cf  a  caloric  equation  tf  state  (27)  and  the  definitions  (30;  and 
(31)  It  is  thi3  relation  n-dch  justifies  phy&icuily  or  use  of  the  words 
nte,r.reratm’sn  and  •  pressure”  for  the  onanfcities  7”  a. ad  &P  as  given  by  the 
fcrml  definitions  (30;  ard  (31).  ‘  * 

In  the  social  cues  when  the  coatinnum  is  hmogeneoue  the  equation  (30) 
become 


<±L 

a  > 

ft 


7l  -  . 

if  1 


(36) 


where  noji  the  diffs rant iai ion  is  along  an  arbitrary  curve  C  in  the  sfs,v 
space.  This  ui%.tonent  is  ouch  stronger  than  tfcf.t  embodied  in  equation  (35) s 
since  the  changes  included  are  net  merely  those  suffered  by  a  given  particle, 
but  my  involve  passage  tc  a  neighboring  p  rticls,,  fortunately  v.e  shall  act 
need  to  use  equation  06)  at.  all,  so  that  <«r  developments  will  be  valid  for 
Inhonogenecua  media* 

Let  uc  now  substitute  the  definition  (19)  of  tbs  tensor  components 
and  'ihe  thcriwyiyasudc  ensrgy  equation  (35)  into  the  aschr.aice.i  crergy  equa¬ 
tion  (26) j 

A*  : -*)/%+£+#%  > 


fi-rjf  tf 


J- 


(37) 


T* 

where  the  digs ipat ion  function  jp  is  given  by  the  definition 
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(33) 


r *  ^ 


•A 


With  the  aid  of  the  equation  of  continuity  (15)  the  equation  (37)  becomes 

fiT§  '(p-r)V*,K  a/%  ,  (39) 

In  the  definition  (19)  of'  the  components  W,  .  of  the  fluid  stress  tenser  the 
scalar  fu.-icti.on  TT  ms  left  arbitrary.  Letruo  now  define  it  to  be  the  thermo- 
dynasde  pressure  p  „  Then  the  energy  equation  (39)  becomes  simply 


/- 


»< 

J(X 


(40) 


Tie  equatiori  (40)  edmits  a  simple  physical  interpretation:  the  rate  at  which 
beat  is  accraulated  by  a  particle  equals  the  rate  at  which  work  is  being  dona 
by  ths  stress  in  deforming  the  continuum  at  the  present  location  of  the 
particle  plus  ths  rate  at  which  heat  is  being  added  to  tho  particle. 


t 


Ths  free  energy  0  is  a  thermodynamic  variable  given  by  the  definition 

(41) 

In  terms  of  0  ths  energy  equation  (26)  beocnen 

fill-  -r  o  {42) 

The  loft  side  is  the  rate  of  increase  of  mechanical  energy/  of  a  given  psrtielo, 
the  right  aids  is  tba  rate  at  which  heat  is  being  token  away  from  that 
particle.  Tc  deduce  the  equation  (42)  we  have  not  used  ths  first  low  of 
themxtynaaics  (27)  or  the  concepts  of  entropy  and  temperature.  The  equation 
(4-2)  in  fact  Is  a  mere  restatement  of  the  mechanical  energy  equation  (26)  t 
dividing  the  internal  onorgr  into  a.  !iaechani.cnl  pn.rt" ,  which  m  may  (Mill  0, 
and  a  "caloric  part,"  which  ve  my  call  sT.  It  is  only  in  the  application 
cf  the  equation  (42)  as  in  pert  IX.,  that  we  use  tho  weaning  of  a  and  T  to 
Snide  no  in  selecting  proper  postulates  for  special  models  of  continue. 

A  continuum  ic  said  tc  be  incog/.jresslble  if  the  density  of  each  particle 
remains  constant: 

p  y  ( % )  4b  )j "  'j  tfo-  J  j  -  (43 ) 

The  caloric  equirc  ton  of  utite  of  such  a  continuum  is  then  of  tho  fora 

'  j  “  r  ./  i  :^f*\)  i  '  *  '  J  J  -  O  (44) 


S8CIA&SIFI3D 
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tte  tSS.  S‘M  in)  18  •*"*»  ««>,  and  .11 

tion  h ts  £:  UUI*  if  p  Is  put  eoual  to  I*  “***«  for  Incosn- 

10n  \-'-5'  becomes  ‘  2ero*  Since  the  continuity  equa- 


V 


<x 

'x 


(4  5} 


3  the  foroe  *“  to  bo  the  ZrlV£ssurly  W°e,,m° 

vr ■  ^SSaUffl-at  nioroojymrio. ,  '  ’ 


irt^^r^rf°var'  r  * « 

“  40  ‘  ™  ^’-^tSg^thot 

(46) 


if 


»/  * 

^  >0<  ~  ^  0«p\ 

for  continuous  madia  in  .  4 

the  Equality  (46)  lot  U  $  le  lad»Pe^ent  of 

',G  in£q,»tlity  (46)  and  the  i^,ra'  *  s 


f/>T§  ^  ^  fm.  , 

eiv.'i'MxK.  s  tj 


(48) 

cmt_an  fDtr°ro  Gretas  but  sith 

8Ur“f  (f  »  »<•  at  leant  o ZmlZ  SLt^SSfT  fr,terlw  “>  "W  <*»„ 
s  -roir  the  surroundings.  Thus  media  .m2  S*  ^  h~at  saJrSy  Hors  into 
*****  reotricttou  oafnrt  ^a“u  tea?  ^*’->»,»»^Utar  &)  holds 
;^'B-1  f"w«r  lot-  aeehanical  ifflf?' '*  ‘*  -cannot  convert 

«taw*at  of  tie  acoond  leu  of  thSvJ!/  ^  that  W)  lo  rot  a  genera’ 
oat;.-  subject  to  the  restriction  f°r  ocntim“'°  >»«»!  It  holds' 

VXI,, 


£S!^lJS3taiS3Ae;hinle..  of  Conte  „■■ 


trary  nertoiuf^„t^hs  fvf"  «*>  behavior  of  an  arbi- 

one  inequality,  '  'reave  thus  obtained  nine  acrtar  aqua  iiena“„d 

Coneei-vatlon  of  fate: 


£4> 


8* 


-h  V 


~  O 


(49) 
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Conservation  cf  Momentum 


rJJ  ~  7~-^ 


(50: 


Conservation  of  energy; 

/$=7 . 


ff\  ft*  ‘  ft  ■ ■  •  ■$?,' e  *.f>)  *0. 


Second  la®  of  Thermodynamics: 


i  - 


(51) 

(52) 

(53) 


(54) 


All  the  above  equations  arc  written  on  the  Eulerian  plan,.  They  form  a  system 
of  9  equations  la  3.8  independent  variables  /O  ,  e,s,  If  as  in 

some  electromagnetic  phenomena  the  force  field  depends  upon  other  fields 
which  are  in  tarn  influenced  by  the  notion  of  the  continuum,  the  number  of 
depsrdent  variables  is  increased,,  So  assume  herb  that  the  components  of  the 
external  field  F-  are  given  functions  of  space  and  time*  Then  to  obtain  a 
determinate  system  we  KUBt  add  9  Bore  equations* 


Ho  further  physical  principles  ore  available,  30  we  must  conclude  that 
the  notion  of  a  vswtonian  continuum  is  in  general  indeterminate.  VJhilo  a 
chosen  specific  foim  for  the  function  f  in  the  equation  of  state  (53)  say 
well  indicate  certain  fundamental  physical  characteristics  of  the  continuum, 
it  cannot  completely  summarise  the  macroscopic  properties  resulting  from 
chemical  composition  and  the  state  of  aggregation  of  tb©  elementary  particles  , 
In  order  to  obtain  a  determinate  system  of  equations  without  resort  to  a  fine 
structure  of  theory  we  gerseralias  our  experience  with  numerous  physical 
bodies  and  construct,  within  the  framework  of  the  mechanics  of  continun  out- 
lln<vd  so  far,  various  p hanomenological  models,  ideal  continun  which  embody 
abstractions  of  ths  behavior  of  classes  of  actual  bodies.  The  simplest  of 
these  models  are  the  elastic  solid  and  the  viscona  fluid,  which  wa  shall 
discuss  presently,  In  defining  a  model  it  13  customary  to  postulate  expres¬ 
sions  for  and  H*  in  tones  of  other  variables.  H1  will  depond  upon  the 
tempers turo  I,  which  dii  not  occur  in  ths  equations  (49)  to  (53),  so  it  will 
bo  necessary  to  -add  the  definition  (30),  Tie  nay  in  faot  eliminate  s  from 
the  caloric  equation  of  state  (27)  with  the  aid  of  the  definition  (30),  then 
©Uri  nate  the  'UJ  jrknes  parameters  as  Infer©,  and  thus  deduce  an  equation  of 
the  form 


Lk  CL  ILs 

SC-J 


(55} 
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iijhich  7,0  may  employ  in  place  of  tho  equation  of  state  (53)  *  The  specification 
of  the  stress  tensor  must  satisfy  the  symmetry  relation  (50)  and  the  inequa¬ 
lity  (54).  We  may  now  eliminate  T1^  and  H*  from  the  equations  (49)  (f>l)  (52) 
and  (55),  which  becomes  a  system  of  six  equations  in  the  six  dependent, 
variables  ,  e,  T,  V*. 

VIII.  The  Classical  Linear  Theories. 


The  classics!  theory  of  elasticity  is  based  upon  a  postulated  linear 
relation  between  stress  and  strain: 


c*«  ?•* 

Jfi  r'<  J 


(56) 


where  the  C>|  are  constants  of  the  material,  while  the  classical  theory  of 
viscous  fluids  is  based  upon  a  postulated  linear  relation  between  viscous 

stress  and  defoliation  rate: 


C 


A .  Cs 

ifi 


(57) 


where  tho  are  constants  of  the  fluid.  There  is  a  complete  formal  ana¬ 
logy  between  uha  two  theories  when  the  theory  of  elasticity  is  simplified  by 
tho  further  assumptions  of  small  displacements,  so  that  the  Eulerian  and 
Lagrangian  description  of  the  deform tion  approximately  coalesce,  and  the 
further  assumption  of  small  displacement  gradients  (or  equivalently,  small 
strain),  so  that  tha  strain  components  may  be  approximated  by  linear  expres¬ 
sion?.  in  the  displacement  gradients.  Then  in  rectangular  Cartesian  coord  !•’ 
nafcss  the  stross-straia  relations  for  isotropic  materials  boccma 


where  u,v,vr,  are  the  cospormts  of  di aria cement  and  X  and  ju  the  elastic 
cool flclentr.  of  Lava.  In  the  classical  theory  of  isotropic  viscous  fluids 
the  vi s  oous-str ess -  de  for  rat  ion  rate  relations  in  rectangular  Cartesian 
(Eulerian  plan)  coordinates  are 


(59) 


where  VJ:,Vy.  are  tJ 
are  the  coeFflcienta  of 


J, 


mpocents  of  the  velocity  rector  and  A  and  /j 

... ..  -  _ Ity.  In  contrast  to  tho  approximate'  equations 

(58) ,  the  equations  (59;  requiro  no  geometric  or  kinenatio  approximations. 


The  cls.orrical  approach  to  social  models  of  contlnua  has  been  strongly 
influenced  by  the  desire  to  obtain  a  final  formulation  consisting  of  simple 
differential  equations  which  in  simple  special  applications  /  be  solved 
explicitly  la  terns  of  simple  and  familiar  6 pedal  functions.  The  availe- 
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bility  of  high-s^ed  computing  machines  and  of  appropriate  techniques  of 
numerical  integration  suggests  that  the  modern  approach  should  emphasize 
correctness  end  generality  rather  than  formal  simplicity.  Postulated  rela¬ 
tions  of  the  form  (58}  and  (60)  obscure  the  conceptual  basis  of  tbs  theories 
of  elasticity  and  fluid  mechanics,  because  they  obtrude  their  linearity  and 
simplicity.  In  this  paper  we  shall  seek  a  more  physical  approach  to  the 
classical  continuum  theories  by  probing  and  clarifying  the  intuitive  physi¬ 
cal  concepts  of  elasticity  and  fluidity  and  then  giving  them  mathematical 
fora.  ?fe  shall  see  that  the  relations  (58)  and  (59)  then  appear,  as  their 
linearity  suggests,  as  first-order  approximations  in  physically  meaningful 
and  inclusive  general  theories. 


A  continuum  in  which  heat  energy  is  not  converted  into  mechanical 
energy  and  mechanics!  energy  is  not  converted  into  heat  energy  will  be  called 
thermally  simple.  The  conservation  of  hoat  energy  is  expressed  by  postulat¬ 
ing  the  equation 

(60) 

It  then  follows  from  the  energy  equation  (42)  that 


In  problems  of  simple  heat  conduction,  defined  by  the  relation  (24),  the 
heat  aquation  (60)  becomes  , 


Only  in  special  cases,  erg»,  the  motion  of  a  perfectly  elastic  body  whose 
elastic  coefficients  are  not  functions  of  T  or  s,  can  the  motion  of  a 
thermally  simple  body  he  determined  without  the  aid  of  equation  (62).  When 
the  motion  is  known,  p  is  known,  and  then  the  equation  (62)  and  the  equa¬ 
tion  of  state 

Ip  it  SlA  &r  „  r~r\  „  («> 

fora  a  system  of  tco  equation©  from  which  the  two  dependent  variables  s  and 
T  may  be  computed,  In  problems  of  heat  conduction  without  motion  the  system 
becomes  simply 


'cUr)  _  At  Y* 

/°  Jr  "  ’•  '*)  j 
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•  •'i.sisCwjw-i 


a! -&I-  „  o  7~)-a  (65) 

P°J  jt* ,  it*,0-'  ■]  it, FJ  jl  J~ 

The  simplest  possible  case  is  that  of  a  homogeneous  body  subject  only 
to  reversible  heat  flow.  Then 

$  -  Sa  •  (66) 

The  equation  (65)  now  becomes  the  ordinary  equation  of  the  conduction  of 
heat: 

'  (67) 

Specifying  a3  a  function  of  the  and  k  as  *  function  of  the  X^  and 
T  then  renders  the  equation  (67)  a  single  diff  re ;,iial  equation  in  a  single 
dependent  variable  T. 

X .  Perfectly  Elastic  Solids 0 

fcany  actual  bodies  bahave  in  a  manner  celled  "perfectly  elastic11  when 
they  are  subjected  to  sufficiently  snail  loads,.  The  body  offers  no  internal 
resistance  to  rigid  body  actions.  The  strain  produced  by  a  given  load  is 
Independent  of  the  rate  at  which  the  load  is  applied,  and  when  the  load  is 
removed,  the  body  returns  to  its  initial  configuration.  Ho  ratter  how  many 
times  loaded  and  unloaded,  its  response  to  the  same  load  is  always  the  same* 
In  short,  an  elastic  body  responds  to  its  present  strain  from  its  preferred 
initial  state  but  exhibits  no  response  whatever  to  its  rate  of  deformation 
or  to  its  pas':  states  of  strain,.  Often  not  stated  explicitly  but  essential 
to  the  definition,  is  the  requirement  that  during  the  deformation  mechanics’ 
and  thermal  energy  are  not  interconvertible. 

fie  shall  non  give  n  precise  mathematical  fora  to  this  somewhat  vague 
concept  of  a  perfectly  elastic  body.  The  energy  requirement  demands  that 
a  perfectly  elastic  body  be  thermally  simple,  so  we  have  equations  (6 0)  nod 
Oil).  The  requirement  that  the  body  respond  only  to  its  strain  from  a  pre¬ 
ferred  initial  state  demands  that  the  free  energy  0  depend  only  upon  tbs 
therrsodynemic  state,  the  components  *a, j  of  the  gradients  of  the  defonaa- 
t:ion,  the  components  j^c  and  g.y  of  the  Lagrangian  and  Euler ian  metric 
tensors  respective?^, 'and  the  ihnomogeneity  parameters  8^: 

f  (*J  ~‘$*j  >/° ,  <cL'J  >  %0j  '  ‘  •'  )  - 

The  occurrence  of  p  gives  a  complete  dependence  upon  the  state,  because  0 
itself  ia  a  state  variable  aud  the  introduction  of  a  second  state  variable 
in  the  functional  "orm  (68)  would  necessarily  reduce  that  relation  to  a 
general  equation  of  state  rather  than  a  defining  assumption  from  else tic 
body,  Th*  temperature  T  would  be  a  acre  natural  representative  of  tho  st-at 
tlvan  the  density  p  ,  but  it  Is  not  apparent  hem  the  subsequent  analysis 
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could  then  be  performed  without  Hurnaghan* r  unnecessarily  restrictive  as¬ 
sumption  that  T  is  held  constant  during  the  motion.  T>&  shall  see  that  we 
nay  employ  Q  now,  yet  by  a  suitable  device  reintroduce  T  later,  when  the 
analysis  ii  complete. 

Since  the  strain  components  £**  are  functions  of  g* *  and  *n,  >  only,  the 
equation  of  continuity  (14)  states  that 

(69) 

The  initial  density  fa  is  a  function  of  the  Lagrangian  variables  *a  only,  so 
it  may  be  regarded  as  one  of  the  parameters  3(^),  Then  by  replacing  ,o  by 
its  functional  form  (69)  in  the  formula  (68)  we  see  that  0  must  have  the 
functional  farm 

4  =  t  &  fijj  >  &U)j  "  j  &(*) )  '  (70) 

Since  $  is  an  Eulerian  scalar  and  the  only  scalar  functions  of  the  tensor 
components  gjj  are  numerical  functions,  the  fora  (70)  must  reduce  to 

f  *  $(<j  c*  j  3(>)  *  $fyy  * fyH))  ,  ^ 

The  energy  equation  (61)  implies  that 

r  0  a//  fit  -O  .  ^ 

In  the  functional  form  (71)  only  the  gradients  xa ,  *  depend  upon  the  time, 
so  that 

a  ~  £(%#) .  (73) 

Sf  st 

From  this  point  Kurnaghan's  analysis  needs  no  modification  in  order  to  re¬ 
main  valid  under  our  slightly  more  general  assumptions.  By  Insertion  of 
the  formula  (73)  into  the  condition  (72)  we  may  show  that  (ref„  a,  pp.246-247) 

*  j  j  Bo)j  Bfr),  * '  'y  B(k)),  (74) 

Hence  again  using  the  energy  equation  (61)  wo  may  finally  deduce  (ref.  a 
pp. 256-257,  247-249)  the  general  Lagrangian  stress-strain  relations! 


iJ  V'4Jfs>  i 


the  general  Eulerian  stress strain  relations: 


where 
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j  ~  <S  %  '«/V  /j 

and  the  stress- strain  relations  for  isotropic  media: 


It  can  be  shown  that  the  last  formula  is  consistent  rrith  the  symmetry  re¬ 
quirements  (17)  (ref.  a,  pp.  248-249).  If  in  the  isotropic  case  the  function 
p0  is  expanded  in  a  power  series  in  the  invariants  (16) , 

ft  i>  -yl,  *J&-nJ;ix  +nls + •  •  ■ • ,  (79) 

where  the  ooefficientso<y\^l,  m,  n, ...»  are  functions  of  the  Lagrangian 
coordinates  only,  the  stress-strain  relations  (73)  become  (ref.  a,  pp.  250- 
25!)  v  - 

-  fix'*  (X-^Tj  ■+  (s£+n\-\’‘Jot  )1,  fa  (Vn-<-  ^ «}T,  Jf/ 

J  *  (80) 

4- [?(#-*) -fa '* °<)£i  J Ej 

4-  a\.Tj  •  •  f 

whore  the  tensor  components  are  computed  from  the  equations 

-  <81) 

In  the  formula  (80)  ne  have  written  down  only  terms  up  to  and  including  those 
of  second  order  in  the  components  Ej  .  If  we  add  the  artificial  recuirr  -> 
ment  that  the  initial  stress  must  vanish,  ire  find  that  cyto  .  Under  the 
hypotheass  of  small  strain  and  small  displacement  the  general  stress-strain 
relations  reduce  to  the  classical  formulas  (58),  which  now  appears  aa  a 
first-order  apjiroxiraation  in  an  exact  theory  which  completely  embodies  the 
notion  of  perfect  elasticity. 

Since  0  is  a  state  variable  there  exists  an  equation  of  state  of  the 


p  z /){<£,  T.  B(,) ,  8(3) /  •  ’  j  B(h})  , 


The  expression  (68)  becomes 


I  J  y°  T>  $0j  ■ '  j  i  $0/ 

So2rv5.ng  this  equation  for  f i  we  find  that 
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9  tCj  cj$uP~j  ‘ j  j  3(<)j  B(>)J  "  3(k)) ,  (84 - 

Thus  0  Is  a  function  of  T.  Sinoa  T  does  not  occur  la  I.,  I2,  or  I?  it 
follows  from  the  expansion  (79)  that  in  general  the  coefficients  ot 
l,m,n  are  functions  of  temperature  as  well  as  of  initial  density  /%  and 
other  inhonogoneities  represented  by  the  B/,\.  The  dependence  of  the  coef¬ 
ficients  of  elasticity  upon  the  temperature  cannot  of  itself  describe  thermo- 
elastic  effects,  which  we  have  expressly  exoluded  by  the  postulate  that  a 
perfectly  elastic  body  is  thermally  simple.  The  effect  of  non-uniform  heat¬ 
ing  upon  a  body  which  cannot  realize  conversion  of  thermal  into  mechanical 
energy  is  rather  to  introduce  an  elastic  inhovogeneity,  so  that  it  is  not 
surprising  that  we  were  able  to  include  these  thermal  effects  in  a  develop¬ 
ment  valid  for  inhomogeneous  bodies. 

By  expressing  the  components  of  the  tensor  Ti^  in  terms  of  other 
variables  we  have  added  6  of  the  9  more  equations  which  we  found  in  part  VII 
to  he  necessary  to  secure  a  determinate  set  of  equations  for  the  motion  of 
a  Newtonian  continuum.  If  the  ooeffloients  of  elasticity  are  not  functions 
of  T  the  equation  of  continuity  (14)  and  the  dynamical  'ki nations  (13)  may 
be  regarded  as  a  systeo^of  four  equations  in  the  four  dependent  variables 
,La,  a(3af//o  or  3r-,x*,.x- ,  ,  according  as  Eulerian  or  Lagrangian  descrip¬ 

tion  is  employed.  The  motion  is  then  determined  from  purely  mechanical 
considerations,  and  once  it  is  known  it  may  be  substituted  into  the  thermal 
energy  equation  (60)  if  we  wish  to  determine  the  temperature  distribution 
as  well,  provided  we  add  the  remaining  3  necessary  equations  by  specifying 
the  heat-flow  field  components  H*  as  functions  of  the  other  variables  and 
employ  a  suitable  equation  of  state  (63).  If  the  ooeffloients  of  elasticity 
depend  on  the  temperature,  however,  the  problem  of  determining  the  motion 
cannot  be  solved  without  incidentally  determining  the  accompanying  tempera¬ 
ture  distribution,  and  all  the  differential  equations  (14)  (IS)  (60)  (63) 
mast  be  solved  simultaneously.  This  division  of  problems  of  elasticity 
into  two  classes  according  as  the  motion  does  or  does  not  depend  upon  the 
thermal  energy  equation  corresponds  exactly  to  Bjerknes*  a  (ref,  f,p,  84) 
division  of  problems  of  fluid  mechanics  into  the  classes  of  "pure  hydrody¬ 
namics"  and  "physical  hydrodynamics", 

XI,  Viscous  Fluids 

In  a  very  rough  way  we  conceive  a  viscous  fluid  as  a  continuum  which 
when  at  rest  experiences  a  state  of  stress  in  which  the  normal  components 
in  all  directions  are  equal  and  in  which  there  ore  no  shearing  components, 
but  which  when  subjected  to  change  of  shape  experiences  "viscous  stresses". 
The  viscous  stresses  are  independent  of  the  history  of  the  fluid.  There 
is  no  preferred  initial  state  to  which  the  fluid  tends  to  return  If  the 
forces  producing  the  motion  be  removed.  It  is  therefore  impossible  that  the 
viscous  stresses  be  functions  of  the  components  of  the  s brain-rate  tensor 
(8),  since  these,  in  turn,  depend  upon  the  components  of  the  strain  tensor. 

In  contrast  to  a  perfectly  elastic  body,  which  responds  only  to  its  strain 
from  a  preferred  initial  state  independently  of  all  its  history  between  that 
state  and  its  present  one,  a  viscous  fluid  responds  only  to  its  present 
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strain  from  the  infinitesimally  preceding  state,  in  other  words  to  its  rate 
of  deformation.  The  rate  of  viscous  dissipation  of  energy  should  theraforo 
be  a  function  only  of  the  components  of  the  deformation  rate  tensor  (?)  and 
possibly  also  of  the  thermodynamic  state,  and  should  vanish  if  the  deforma¬ 
tion  rate  components  vanish.  i7e  nay  pat  the  preceding  physical  considers- 
tions  into  a  precise  nsathematical  form  thus: 

£ '  $  (e/j  ^0)i  "y  (S5) 

ft/tj-O  if  cJHt*  £j=0 .  (36) 

The  dissipation  function  <&  is  not  a  state  variable,  so  two  state  variables 
must  occur  in  it  in  order  that  its  dependence  upon  the  thermodynamic  state 
be  as  general  as  possible,  rfe  have  chosen  the  temperature  T  and  the  thermo¬ 
dynamic  pressure  as  a  matter  of  convenience.  As  we  stated. earlier,  the 
function  7T  in  the  definition  (19)  of  the  tensor  components  is  to  be  taken 
&b  the  thermodynamic  pressure  p  for  compressible  fluids,  and  for  incompress¬ 
ible  fluids  left  undefined,  Ha  have  included  the  Btress  pressure  p  in  the 
form  (85)  so  that  even  in  the  incompressible  case,  whan  =0,  some  sort  of 
pressure  will  appear  in  0 

Viscosity  m3  discussed  by  Maxwell  in  terms  of  the  following  ideal  ex¬ 
periment.  Let  a  viscous  fluid  be-  confined  within  a  channel  bounded  by  two 
parallel  infinite  walls,  one  of  which  is  moving  relatively  to  the  other  with 
a  constant  velocity.  The  fluid  will  adhere  to  the  walls  and  the  retarding 
force  will  be  proportional  to  the  relative  velocity  and  to  the  ares  in  con¬ 
tact,  and  inversely  proportional  to  the  distance  between  the  walls.  That  is, 
there  exists  a  modulus  such  that 

Fores  on  11  -  A  («•  *  contact).!*) 

• 

We  shall  not  employ  the  relation  (87)  &3  an  actual  algebraic  formula,  nor 
shall  we  utilize  the  ideal  experiment  just  described  to  assist  us  in  our  def¬ 
inition  of  a  viscous  fluid  in  any  way  except  to  postulate: 

There  exists  a  material  constant  yU0  of  the  dimensions 

C/J  8 

The  postulates  (85)  (86)  and  (88)  we  shall  take  as  the  defining  postu¬ 
lates  of  a  viscous  fluid. 

Let  <£  be  expanded  in  a  power  series: 
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fix.  *  '  • 


where  the  coefficients  ara  of  the  forts 


•  ■  .  •  ,  .  (S?) 

«a  *  j*iab  ,  %  *>)• 


- 1,  <v  ••(«-» 

./•jv*  •)**-» 

The  postulate  (36)  and  ths  formula  (33)  then  imply  that 


(91) 


For  the  succeeding  dimensional  argument,  let  us  suppose  all  formulas 
taken  in  rectangular  Cartesian  coordinates,  so  that  all  tensor  components 
assume  their  proper  physical  dimensions.  Then  since 


DO*  jpfey-jtisrp  ’ 

it  follows  from  formula  (89)  that 


(92) 


i, t%gg-3 _ Ljanklfl! 

length  J 


(93) 


From  formula  (90) 

c"  "•  fr)) 


(94) 


Then  by  the  postulate  (88) 


*  * 

LJ z(  ( b , pJ Tj B(,x B(k) )j 

JAt  » 


(95) 


where  the  function  -f  is  dimensionless.  Th®  parameters  B/jq  we  shall  take 
dimans ionless,  &3  if  they  ara  not  dimensionless  they  may  be  divided  by  ft 
reference  quantity  of  the  same  dimension.  Than  since  the  dimensionless 
relation  (95)  connects  three  quantities  of  two  indapandent  dimensions,  it 
must  not  involve  the  temperature  T  at  all.  unless  there  exists  a  reference 
temperature  77  and  then  it  must;  be  of  the  form 
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Sow  ‘Yotr  fortula  (90) 


<1&  «i,h  fa  TiA>>  &>«  9(W'  - 


(97) 


Sot  counting  quantities  that  are  diEcnsionless,  we  nay  regard  this  formula 
as  a  relation  among  the  3^  quantities  Cjy  of  dire  as  ion  [base]]/  [length] 
[tiraej,  the  quantity  of  dissension  finasu  J/  [length J,  and  the 

pressure  p  of  dimension  \js&ss  ]  /  [length  j[tirej z .  Since  theBe  Tr  2 
quantities  are  all  exprasaible  in  terns  of  the  two  fundarental  dimensions 
i_n» 3'!] /I. e eg th  J  »>>d  ftinejf.  tfctre  rust  be  a  relation  among  3r  dimension¬ 
less  ration  formed  from  then: 

(Ml ) .  0  m 

MtUJ  ' 

8y  the  f  crania  (96)  rro  nay  express  this  result  ir.  the  form 


f 


GtL  ~7T"^(  J  )~f6 } 


(99) 


where  to  have  omitted  to  rrite  the  3^  quantities  i Ufa  explicitly  in  the 
function  {  at  the  right  since  they  are  dimensionless  functions  of  the  other 
arguments  of  £  „  Thi3  reasoning  can  be  repeated  for  the  coefficients  of 
each  order,  so  that  we  obtain  aa  our  final  fora  for  the  dissipation  function 


<«■/ 


r  <.  ^.3—  ifUM-s-vm  -A 

$  f"’*  /3>A'"A  ^ 


(100) 


where 


= '’v-a  •••/»«  A  a«) 


and  is  diE3ncionlo3p, 

A  fluid  is  isotropic  if  $  is  invariant  in  form  under  transformations 
of  ths  Eulerian  coordinates,,  Since  any  invariant  function  of  a  3  rowad  square 
matrix  is  a  function  of  three  independent  invariants  of  that  ratrix,  and  since 
thr  quantities:  ,  _  c 

form  such  a  3ot  of  invariants,  for  isotropic  fluids  £  pust  be  of  tlie  form: 

■*■/$(*&*$*  e~* 

(102) 


UNCLASSIFIED 


-22° 


NCLM  9223 


where  /  Kp  Kg,  are  dimensionless  functions  of 

p>  /p,  T,4T0,  s(2)j,**»  B^\.  >  and we  shall  call  "first-order  co¬ 

efficients  of  viscosity,"  Ki,  Kg,  K~  "second-order  coefficients  of  vis¬ 
cosity",  etc.  Of  course  the  coefficients  of  viscosity  A  and  yU  used  here 
have  no  relation  to  the  constants  of  Lame  #A  and  ju  used  in  the  theory  of 
elasticity. 

The  expansion  (102)  is  tho  analogue  for  viscous  fluids  of  Viurnaghan's 
formula  (79)  for  perfectly  elastic  solids.  With  the  aid  of  the  definition 
(33)  we  may  deduce  from  it  the  following  viscous-stress  -  deformation-rate 
relations ;  a 


C  -X  £*  £j  +■  &■  +  ‘j-  £|C(C/<f/ 

^1:1; +  K3£i  £*)+  ■■■  ■ 


(103) 


The  classical  Kavier-Stokes  equation  (59)  thus  appears  as  a  first-order 
linear  approximation.  In  the  theory  of  elasticity  the  criterion  for  break¬ 
ing  off  the  expansion  (79)  after  the  quadratic  terms  is  simply  that  the 
strain  be  small.  To  formulate  an  analogous  criterion  for  viscous  fluids 
let  ua  introduce  the  dimensionless  number  % \ 


m  *s  . 


<104.) 


Then  the  linear  llavier-Stokes  equation  is  a  sufficient  approximation  if 


lr  <  <.  x 


(105) 


Thus  the  classical  linear  theory  may  bo  expected  to  become  inadequate  for 
very  viscous  fluids,  for  fluids  subjected  to  very  low  pressures,  or  for 
fluids  being  very  violently  deformed.  If  we  suppose  that  X  ,  then  for 
water  at  a  stress-pressure  of  one  atmosphere  the  enormous  distortion  rate 
of  2,000,000  sec'1  is  required  in  order  that  IT  bo  of  the  order  of  1. 
Burnett's  derivation  of  the  equations  of  slip  flow  from  kinetic  theory  gives 
vi3cous«otreso  -  deformation- rate  relations  which  essentially  agree  with  the 
formula  (103),  as  far  as  the  latter  goes,  but  Include  also  terms  involving 
gradients  of  the  thermodynamic  state  which  the  present  theory  does  not  in¬ 
clude,  since  they  do  not  satisfy  the  defining  postulate  (85).  Burnett's 
result  indicates  that  when  a  gas  becomes  sufficiently  rarefied  that  it 
ceases  to  obey  the  Havier-Stokes  equations  it  ceases  also  to  exhibit  some 
of  the  fundamental  properties  of  fluids  (ref.  b). 

From  the  formula  (103)  we  see  that 

<  .  fa  +•  +  f\(3K,  f  K,)i,  \  K,  ?;  if]+  ■ .  .  (106) 
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Then  from  the  definitions  (19)  and  (22)  since  7T-©  it  follows  that  for 
compressible  fluids 

3(P -/>)-.  (1M iu)f,  - Aj  E;  £?><•  -(107) 

Thus  the  thermodynamic  pressure  P  and  the  stress  pressure  p  are  in 
general  unequal 0  In  order  that  the  coefficients  of  each  term  in  the 
expansion  (107)  nu3t  vanish:  r  ’ 


3\  +1/J  -O  j 

3k ,  +  K x  —Oj 

Kj  =  <V 


(108) 


The  first  of  those  conditions  on  the  material  constants  of  the  fluid  ie  the 
classical  Stokes  relation.  It  easy  to  see  that  in  general  the  number  of 
independent  reorder  coefficients  of  viscosity  ia  equal  to  the  number  of 
different  partitions  of  r  as  a  sum  of  multiples  of  1,  2,  and  3,  but  if 
p-p  there  will  remain  only  one  independent  coefficient  of  viscosity  of 
each  order.  Dubes  (quoted  in  ref,  g,  p.  498}  has  shown  that  for  the 
classical  linear  theory  the  second  law  of  thermodynamics  (48)  is  equivalent  to  t 
to  the  conditions 


3\+-2/t£C 


(109) 


His  analysis  is  not  easily  extended  to  the  general  case,  and  it  is  not 
dear  what  are  the  proper  generalizations  of  these  inequalities.  For  in¬ 
compressible  fluids  the  function  TT  has  not  been  specified  and  the  continu¬ 
ity  equation  (45)  implies  that  Xt~0  ♦  Instead  of  equation  (107)  we  have 


X-p)  --  ^r 


K. 


£/S  oi  4 


(no) 


Hence  in  the  classical  linear  theory  of  incompressible  fluids  even  though 
the  function  7T  in  the  definition  (19)  was  left  porfeotly  general  the 
nature  of  the  Navier»Stokes  equation  forces  it  to  be  the  stress  pressure, 
but  in  the  general  theory  the  equality  TT  -  p  is  equivalent  to  a  sequence 
of  relations  among  the  coefficients: 


an) 


XII,  Conclusion 


The  familiar  analogy  between  the  theory  of  elasticity  and  tbs  theory 
of  viscous  fluids  rests  upon  the  rather  superficial  similarity  of  the  class¬ 
ical  linear  force  laws  (58)  and  (59)*  A  more  physical  and  less  heuristic 
formulation  of  these  two  models  of  Newtonian  continue  reveals  them  to  be 
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rather  opposites  than  analogues,,  kn  elastic  body  is  a  conservative  mechan¬ 
ical  system  which  responds  only  to  present  strain  from  a  preferred  initial 
state  and  which  cannot  interconvert  mechanical  and  thermal  energy,  while  a 
viscous  fluid  is  a  dissipative  mechanical  system  in  which  the  rate  of  dis¬ 
sipation  of  energy  depends  only  upon  the  instantaneous  rate  of  deformation. 
Khan  these  concepts  are  expressed  in  mathematical  terms  we  obtain  the 
general  non-linear  fores  laws  (30)  and  (102),  respectively,  between  which 
no  formal  analogy  is  apparent.  The  two  theories  now  exhibit  a  deeper  kin¬ 
ship,  however,  in  that  each  is  the  tneory  of  a  special  model  of  a  Newtonian 
continuum  subject  to  all  the  general  laws  listed  in  part  VII,  and  each  is 
defined  by  specifying  the  functional  dependence  upon  suitable  variables  of  a 
scalar  function  representing  tha  energy  or  rate  of  dissipation  of  energy,, 
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